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Abstract 

We derive the AC Josephson frequency using the gauge invariant equations of motion with 
including the battery-contact effect. The frequency is given as a sum of a contribution from 
an Aharonov-Bohm phase that arises when charged-partices pass through an electric field in the 
insulator and that from a chemical potential difference arising form the battery contact; each of 
them contributes |g|V//i to the frequency where q is the charge, thus, the sum is 2\q\V/h. The 
observed Josepshon frequency, 2\e\V/h, hence, means that the charge on the tunneling particles is 
e. A variety of derivations so far miss one of the above two contributions; the original derivation 
misses the first contribution due to the lack of inclusion of the electric field in the insulator. The 
present result indicates that the phase of the order-parameter does not arise from the number 
fluctuations of Cooper pairs. We present an alternative origin for it; we argue that a Berry phase 
that arises from spin-vortices is the origin. 

PACS numbers: 74.50. +r,03.65.Vf,74.20.-z,74.72.-h 
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I. INTRODUCTION 



The origin of persistent current in superconductors is believed to be a solved problem^. 
One of the foundations of the current understanding of the persistent current generation 
is the appearance of an AC current of a frequency 2\e\V/h in a superconductor-insulator- 
superconductor junction with an application of a constant voltage V, where e is the electron 
charge and h Planck's constant. The prediction of this effect (the AC Josephson effect^) 
and its observation 3 established the validity of the current understanding. Especially, 2|e| 
in the Josephson frequency 2\e\V/h has been considered as the direct evidence that paired- 
electrons tunnel through the junction in accordance with the BCS theory^. In this work, 
however, we show that the frequency 2\e\V/h actually indicates that the persistent current 
carriers are not paired-electrons but electrons. 

A variety of derivations so far have defects. A common defect originates from the usage 
of the Cohen-Falicov-Philips tunneling Hamiltonian^ without including the electromagnetic 
field in the insulator.- In the Josephson tunneling case, the voltage-drop occurs across the 
insulator between the two superconductors of the junction, and an electric field that bal- 
ances the battery voltage exists in the insulator; when a charged particle passes through 
the insulator, it is either accelerated or decelerated by the electric field. This situation is 
markedly different from the ordinary tunneling case where the insulator acts simply as a 
potential barrier. 

The preset result calls for rethinking of the origin of the persistent current generation 
in superconductors; it indicates that Cooper pairs play only a secondary role as to the 
persistent current generation is concerned. Since the bulk properties of superconductors are 
well-explained by the Cooper pair formation, an explanation that is compatible with the 
above observation is that tunneling particles are broken-Cooper-pairs or unpaired electrons 
that exist in the surface region. In this respect, it is worth noting that serious deviations from 
the BCS theory have been observed in the Knight shift measurements; and such deviations 
have been thought to reflect the difference of electronic states in the bulk and in the surface 
region^. 

As an alternative for the conventional persistent current generation mechanism, we will 
put forward a novel one that has been proposed, recently 9,10 . It utilizes a Berry phase 
arising from spin-vortices^; the persistent current is realized as a coherent collection of 
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stable loop currents that are generated around spin-vortices by the fictitious magnetic field 
of the Berry phase.- In the latter half of the present work, we will show that this mechanism 
(we will call it the "spin- vortex superconductivity") explains the observed AC Josephson 
frequency 2\e\V/hr-, the flux quantization in the unit hc/2\e\ (c is the speed of light and 
the oscillation pattern of the SQUID 13 . 

If the spin-vortex superconductivity is the right one, the origin of the phase of the or- 
der parameter is the Berry phase arising from spin- vortices. In this respect, it is worth 
noting that there are a number of indications that spin-vortices are present in cuprate 
superconductors^; thus, the present view may provide a new perspective that leads to the 
elucidation of the mechanism of the cuprate superconductivity. 

II. A GAUGE INVARIANT DERIVATION OF THE AC JOSEPHSON FRE- 
QUENCY 

Let us derive the AC Josephson frequency in a gauge invariant manner using equations 
of motion for the order parameter. 

First, we denote the superconducting order parameter as 

vl/ = p^e", (1) 

where p is the supercurrent carrier density and 9 is the phase of the order parameter. In 
the conventional theory p is the number density of Cooper pairs, and 9 is the phase variable 
that is canonical conjugate to p. 

We consider an SIS (superconductor-insulator-superconductor) junction which is con- 
nected to a battery of the voltage V as seen in Fig. []Ja). In an ideal situation, the voltage- 
drop that balances the battery voltage occurs entirely across the insulator between the two 
superconductors, and the resistance is zero. This means that the junction is a capacitor in 
the zeroth order approximation, thus, surface charges develop so that the generated electric 
field balances the battery voltage. 

The Josephson current is given by 

J = Jq sin 0, (2) 
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where Jo is a constant, and is expressed as 

*=/ , (s A - v 0-*' (3) 

h is Planck's constant divided by 2tt, q is the charge on the current carrier, and A is the 
vector potential; the integration is performed along a path that connects points in the 
superconductors Si and S2 through the insulator I. 
The combination 

W-£A (4) 
is gauge invariant since for the gauge transformation 

A' = A + V/, (5) 

9 transforms as 

* = » + (6) 

where / is a single- valued function. The current in Eq. (j2J) is gauge invariant, accordingly. 
The following combination is also gauge invariant 

*+f (7) 

because the scalar potential (p transforms as tp'=(p—f/c. In the following, the gauge invariant 
combination in Eq. (J7|) plays a crucial role. 

The basic equations of superffuid phenomena^ are given by 

1 sn 

where TC is an energy functional. 

For the superconductivity case, the above equations should be modified as 

p = nW (10) 

■ q 1 5Ti. 

9+ ? = iv < n ' 
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FIG. 1: (a) A charged particle with charge q tunneling through an superconductor- insulator- 
superconductor junction that is connected to a battery of the voltage V. The junction is a capacitor 
in the zeroth order approximation, thus, surface charges develop so that the resulting electric field 
balance the voltage V. The particle tunneling from Si to S2 is accompanied by another charge- 
transfer from S2 to Si through the battery which fills the vacant level created by the tunneling 
electron, (b) A charged particle with charge q passes through the capacitor of the voltage V; it 
acquires the phase shift corresponds to the energy gain of \q\V. 



to take into account the interaction of electromagnetic field and charged-particles. Here, 
H is a gauge invariant energy functional, H = H[p, V6> — ^A], and ^ is a gauge invariant 
chemical potential, in which V(9 appears in the gauge invariant combination V6* — ^-A. 

Let us calculate the time-derivative of 0. Using Eqs. (j3J) and (TTTI) . and the relation 
E=— d t A/c— V</?, we obtain 

r-2 







m 




2 





(12) 



There are two contributions to it. The first is a contribution that arises when a charged- 
particle passes through the electric field in the insulator (Fig. [TJ(b)). Since the electric field 
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by the surface charges balances the voltage of the battery, it is calculated as 

< 13 » 

This is the Aharonov-Bohm phase contribution. 15 This term is absent in the original 
derivation;-^ there, the insulator is taken as a potential barrier just like an ordinary tunnel- 
ing problem.- 

The second is a contribution from the gauge invariant chemical potential difference pro- 
duced by the battery contact. Since the battery maintains the chemical potential difference 
qV, it is calculated as 

1 (5H 5H 

h\5p 2 5p 

The constant voltage is maintained through the charge flow into and out of the junction (see 
Fig. 11(a)); thereby, the work done by the battery is added to the electronic system in the 
junction. This contribution will be missed if the junction is treated as an isolated system. 

Overall, the time-derivative of <ft is given by <fi = 2qV/h, which corresponds to the Joseph- 
son frequency 

The observed frequency^ is 2\e\V/h, thus, q should be e. This means that the charge on the 
tunneling particle is e, not 2e. This reverses the current understanding that the Cooper-pair 
tunneling is responsible for /j. 

In order to clarify the above rather abstract argument, let us examine the same problem 
from the energy point of view by identifying the charged particle as an electron. In Fig. [21 
the energetics for the tunneling process is depicted. The chemical potentials in Si and S2 
are and /i 2 , respectively. When the tunneling process is absent, the battery creates the 
situation where all states below fi\ and /i 2 are occupied in Si and in S 2 , respectively. 

When the tunneling from Si to S 2 through the insulator occurs, following events happen, 
roughly: First, the tunneling-electron is accelerated by the electric field inside the insulator, 
and gains the energy \e\V (this energy is taken from the energy stored in the electromagnetic 
filed in the insulator). Second, the hole in Si created by the tunneling is filled by an electron 
that is transferred from S 2 to Si through the battery; the energy gain by this process is 
A*i ~ A*2 — which is the work done by the battery. The total energy gain is, thus, 




FIG. 2: Energy diagram for an electron tunneling through the SIS junction in Fig. [TJ hi and fi2 are 
chemical potentials in Si and 52, respectively. One of the most simple tunneling cases is depicted; 
initially, the tunneling-electron has the energy jXi- It gains the energy \e\V from the acceleration 
by the electric field in I. The electron tunneling is accompanied by an electron transfer from S2 
to Si via the battery; this gives rise to the energy gain Hi — pb2- The total energy gain for the 
electronic system by this process is, thus, \e\V + fii — [i% = 2|e|V. When this energy is discarded, 
the transfered electron fills the hole in S2 created by the electron transfer through the battery. 
Black and red electrons schematically describe electrons before and after the process, respectively. 

\e\V + fii — fi2 = 2\e\V. Third, by deexcitation of the energy 2|e|V, the tunneled-electron 
fills a hole in S2 that is created by the via-battery-electron-transfer. If the discarded energy 
is used for an excitation of an electron in S2, a reverse tunneling process will start. Actually, 
there are many ways to fill the holes created by the tunneling and the via-battery-electron- 
transfer; thus, the real processes are much more complicated; nevertheless, the net energy 
gain or loss is 2\e\V for each tunneling process. 



III. MODELS FOR THE JOSEPHSON JUNCTION 

The result obtained in the previous section is surprising; it reverses the current under- 
standing that the Cooper-pair-tunneling is the origin of the observed AC Josephson fre- 
quency. Here, we consider it, again, using a textbook model 16 . This will make the general 
argument given above clearer. 
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The system of equations of the model is given by 

iTiJ^tfi = T^ + zn^ (16) 

ih^* 2 = + (17) 

where = pj^ 2 e* 6 ' 1 and \&2 = p^e* 6 * 2 are order-parameters in Si and S 2 , respectively; T is a 
transfer integral between the two superconductors through the insulator; and p,\ and p 2 are 
chemical potentials of Si and S2, respectively. The battery produces the chemical potential 
difference /i 2 — p\ = qV; this gives rise to the contribution that corresponds to Eq. (TTJJ. 

When a battery is connected to the junction, the equations in Eqs. ffTB"]) and ffTTl) have 
to be modified; in order to include the electric field generated by the surface charges, scalar 
and vector potentials must be introduced, 

d 2 

ih — #1 = Te*£ & A dr ^ 2 + q^i^i + Pi^i (18) 
ot 

ihft*2 = Te 1 ^ £ Adr ^i + g<p 2 * 2 + /i 2 ^ 2 . (19) 

Actually, the inclusion of the vector and scalar potentials is crucial to obtain gauge invariant 
results, as seen below. 

It should be remembered that Eqs. ( Tl8|) and ( fl9l do not include the charge flow into and 
out of the junction; such a process has to be included "by hand". 

From Eqs. ([TBI) and (fl9l . equations for pi,p 2 ,9i, and 9 2 are obtained; 



Pi 



-P2 = ^(pip 2 ) 1/2 sin (j\±A-V0)-d^ (20) 



rp / \ 1/2 / z-2 

9 Pi T p 2 \ 1 ( f . q 



m / \ 1/2 / /.2 

g p 2 T px\ 1 f f q 



* + 5* = -T-ft £ o / (^A-Vfl.*). (22) 



The r.h.s. of Eq. (120]) corresponds to Eq. (jzj); actually, if we take into account the charge 
flow into and out of the junction "by hand" as depicted in Fig. [HI the condition pi=p 2 =0 
is established. The second and third correspond to Eq. (jlip ; by imposing the condition 
Pi = P2 "by hand", the gauge invariant chemical potential difference across the junction is 
maintained to be qV. Then, the relation in Eq. (fl4j) is obtained. Thus, the final result is 
Eq. (US). 



8 



Sj I s 2 




FIG. 3: The charge flow between the two superconductors (solid-line arrow), and those added "by 
hand" (dotted- line arrows) to take into account the battery contact effect. 




(a) (b) 



FIG. 4: Schematic plots of density of states, (a) the SIS junction biased by a voltage V that allows 
tunneling from the filled band in Si to the empty band in S 2 - The filled and empty bands are 
colored in blue and yellow, respectively. The DC current flow with dissipation occurs. The junction 
is like a resistor in this case; (b) the same as in (a) but the voltage is not enough to allow tunneling 
from the filled band in Si to the empty band in S2. The voltage of the battery is supported by the 
surface charges developed in the interfaces. The junction is like a capacitor in this case. 

It may be worth noting that /j cannot be replaced by g<£>; /J is gauge invariant but qip is 
not. The gauge invariant chemical potential difference brought about by the battery-contact 
has to be included using //; this point becomes apparent if we adopt the gauge in which ip = 
and take the limit T — > in Eqs. (fl~8l) and (fT9l) ; then, the only way to include the chemical 
potential difference is to use //i and fi 2 - 

A very often used model for the SIS junction is the Cohen-Falicov-Phillips tunneling 
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Hamiltonian- given by 



if C FP = H L + H R + H T 



(23) 



where Hi,, Hr, and Ht are the Hamiltonians for the left superconductor Si, right supercon- 
ductor S2, and the tunneling through the insulator, respectively. 

Two types of tunneling currents in the SIS junction are currently explained the above 
Hamiltonian. One is the single-particle (or quasiparticle) tunneling current, and the other 
is the current generated by the collective mode 9 (usually, it is called, " the Cooper-pair 
tunneling", but we will not use it, since it is reversed in the preset work). Both of them are 
usually calculated by taking as a perturbation. 

The single particle current is dominant in the situation given in Fig. H^a). This yields the 
current that reflects the density of states of the two superconductors including the chemical 
potential difference brought about by the battery-contanct;— 



where I and r denotes the quasiparticle states in Si and S2, respectively; their energies are 
Ei and E r , respectively. T\ r are tunneling matrix elements, and // and f r are distribution 
functions for the particles in the states I and r, respectively. 

Here, the effect that corresponds to Eq. (fl4j) is taken into account by the energy shift eV 
in the argument of the delta function. In this case, a DC current flows, and the junction is 
more like a resistor than a capacitor. The phase 9 is not definable throughout the junction. 
The matrix element T/ r simply describes the potential barrier penetration rate through the 
insulator; thus, the phase effect that corresponds to Eq. (fl3l) does not arise. 

On the other hand, the other type of current will be dominant in the situation given in 
Fig. Htb), in which the phase 9 is definable throughout the junction; this corresponds to the 
Josephson effect case. The junction is like a capacitor and an AC current flows (the DC 
Josephson case should be viewed as the V — > limit of the AC case). In the conventional 
theory, 9 arises as the angle-variable that is canonical conjugate to the Cooper-pair number 
density. The angle-variable is included in the BCS wave function as 




(24) 



l,r 




(25) 



k 
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where the parameters satisfy + v\ = 1. 9 has been thought to be canonical conjugate 
to the number density of Cooper-pairs, thus, the Josephson current is often said to be the 
Cooper-pair tunneling current. 

Temporal and spatial variations of 9 with added vector and scalar potentials give rise to 
the additional contribution given in Eq. (IT5|) : the magnetic field counterpart caused by the 
spatial variation of 9 with the added vector potential is known to give rise to a very sensitive 
current flow to an applied magnetic field as is seen in the SQUID experiment. 

In the original derivation for the AC Josephson frequency, the contribution given in 
Eq. ( Tl3l) is missing;- 1 ^ this may have occurred by the direct use of the theory that is suitable 
for the quasiparticle tunneling case to the Josephson tunneling case;- then, the phase factor 
for the Aharonov-Bohm phase was not included. There are other types of defects in the 
derivations: sometimes, the substitution of \x (which is gauge invariant) by qcf) (which is 
gauge covariant) occurs, tacitly; 18 in this case, the contribution that corresponds to Eq. (fl4|) 
disappears. Yet, another mistake is related to treating the junction as an isolated system; 
then, the contribution that corresponds to Eq. (TH1) also disappears.— 

IV. SPIN- VORTEX SUPERCONDUCTIVITY 

In the conventional superconductivity theory, the phase of the order parameter, 9, is 
thought to arise from the fluctuation of the number of Cooper-pairs^ 4 . In this interpretation, 
however, the charge on the tunneling particle is 2e, thus, contradicts the experimental result 
according to arguments in previous sections. In the following, we present an alternative 
origin of the order parameter. Actually, we will show that an electronic system with spin- 
vortices has 9 and p that satisfy the equations of motion in Eqs. fflUj) and ffTTJ) with the 
e-charge carrier. 

Let us derive the phase 9 that appears in the system with spin-vortices. 
We consider the Hamiltonian given by 

where the first term describes electron hopping, the second term is that of the scalar poten- 
tial, and the third term is the rest. 
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When spin-vortices are present the following new electron operators cij and bj are 
convenient;— they are related to the original operators, and Cji as 



hi V2 \ -e 



i ' - \ —e 2 e 2 




(27) 



Spin operators at the jth site are given by 



1 i 
S x (j) = - cos ^(a]aj -b]bj) + -sin ^(a]bj -b]aj), (28) 

1 ■ i 

= (3°) 
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If we express expectation values of a^aj — b^bj and z(ajfej — &]%) as (a^aj — b^bj) =Aj and 
Ojbj — b^aj)—Bj, respectively, we have 



<^(j)> = jAf + Sfcosfe-e*,), (31) 



(5 y (;)> = y.l] • Bjsinfo- aj). (32) 

where «j = tan -1 Bj/Aj. Above equations indicate that the phase £j—otj is the spin direction 
in the x-y plane, thus, £ and £ + 2n are physically equivalent. The phase factor exp(ixj/2) 
in Eq. (127]) is introduced to ensure the single-valuedness of the transformation matrix by 
compensating the sign-change of exp(±z£j/2) that occurs for the shift £ — > £ + 2n. We may 
take x — 6> but other choices are also possible. 

For the gauge transformation in Eq. (jSJ), electron operators are modified as 

c ja -»• exp(-ij-//) c i(T , (33) 

which, according to Eq. ( 1271) . means that x is modified as 

xj = - (34) 

therefore we may take 9 in Eq. (ED to be 9 — — x/2. 

Using the new creation and annihilation operators, the hopping term in the Hamiltonian 
Eq. (126|) is now written as 

-XV^ A4 ^ Hr [cos ^(4a i+ 6t 6j) + i S in ^(afe+6^] . (35) 
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Note that A and Vx appear in the gauge invariant combination (qA/hc + Vx/2). 

Now, let us derive equations of motion for the phase \ and the number density p by 
employing the time-dependent variational principle with the following Lagrangian, 

£ = <*| - ffW (36) 



where, the state vector is expressed as 

i*> = X;/«i«>; ( 3? ) 

as 

ot — ( n ai, ■ ■ ■ , n aNsi n bi, 4 ' ' , ^6Ar s ) denotes an occupation pattern of electrons, where N s is 
the number of sites, n a j and n b j are eigenvalues of the number operators a]aj and bjbj, 
respectively. The state vector \a) is defined as 

i«) = n a ) n 6 ii yac )' ( 3g ) 

jess kesg 

where S% and S£ denote sets of sites where the occupation numbers in ajcij and b^bj are 1, 
respectively. 

Noting that a] and fe] are time-dependent because Xj an d Cj are time-depedent, the 
equation (1361) yields, 

C=ihJ2fah+J2^^-J2f*M^\H\P)^ (39) 

a j ot,f3 

where terms with the time derivative of £ cancel out. 

From the above Lagrangian, the Hamiltonian is obtained as 

n=Y,rJMH\P)- (40) 

o,/3 

Now, consider the ground state for given V£, V%, A, and tp. The Hohenberg-Kohn 
theorem^ tells us that it is a functional of p; then, we may express it as 

a,P j 

where the Coulomb potential term is separately treated for the later convenience. 

Next, we make a crucial assumption: we consider the situation where low energy excita- 
tions that give rise to the time dependence of f a are absent; this means that the interaction 
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Hamiltonian H mt is such that it gives rise to an energy gap that suppresses the time de- 
pendence of f a . Actually, this is indeed what is observed in the model considered in our 
previous work 10 . The energy gap formation by the BCS theory may be considered as such 
a kind. 

Then, low energy excitations are described by the following Lagrangian: 



Tic 

p,A+— Vx, V£ 

2( 1 



The momentum conjugate to Xj is, thus, obtained as 

5C _ dC % 



this shows that p and — x/2 are conjugate variables. 
Finally, the equations of motion are given by 



Pf 



P 



_26H 

h 8x 

2 5H 



. _ 2g = 

h h 5p 



(42) 



(43) 



(44) 
(45) 



which correspond to Eqs. (|10p and (fTTl) with 9 = — x/2. 



V. FLUX QUANTIZATION IN THE SPIN- VORTEX SUPERCONDUCTIVITY 

It has been thought that the quantization of the flux in the unit $o — hc/2\e\^- is the 
evidence that the persistent current carriers are Cooper pairs. However, we will show that 
the same flux quantization is explained by the spin- vortex superconductivity. 

In the spin-vortex superconductivity, only V% + ||A and p are active variables in low 
energy excitation. In the ground state, it is expected that Vx + f| A is zero and the electric 
current is zero. If the Meissner effect occurs, the electric current is absent deep inside the 
sample; then, the condition Vx + I 2 A = should be established deep inside the sample. 

Let us consider a ring-shaped sample and consider the line-integration of A along a path 
that is deep inside the sample. If this path is so chosen that it encircles the hole of the ring, 
we have 

j) A ■ dr = j> V% • dr = ^p| n (46) 
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since x is an angular variable with period 2ir; here, q = — |e| is substituted. The above for- 
mula shows that the flux is quantized in the unit $o — hc/2\e\ as is observed in experiments.— 
The assumed Meissner effect above will arise from the stability against the fluctuation 
of Vx + A; a nonzero Vx + A should increase the total energy, thus, expelling of 
the magnetic field from the bulk by keeping Vx + |fA = is energetically favorable. This 
expelling of the magnetic field is nothing but the Meissner effect. 



VI. AC JOSEPHSON EFFECT IN THE SPIN- VORTEX SUPERCONDUCTIVITY 

Now we consider the AC Josephson effect in the spin-vortex superconductivity. 
By introducing the "order parameter" ty s defined as 

the Lagrangain in Eq. fflSl) is expressed as 

C = [ d 3 r^* s (ih^- -h-qip] ^ s (48) 



dt 

where h is defined as the operator that satisfies 



Tic 

p,A+— Vx,V£ 
2q 



(49) 



and the conservation of the number of particles J d 3 rp = is used. 

From the stationary condition of the Lagrangian in Eq. (T4*8|) with respect to the variation 
of we obtain a "Schrodinger equation", 

in—m s = (h + qip]m s (50) 



dt 

where is the "wave function", and (h + qip) is the "Hamiltonian" . 

For the Josephson junction problem, we introduce supports |1) and |2) in the supercon- 
ductors S\ and 5*2, respectively; they satisfy 

1 = (1|1) = (2|2), 

= (1|2) = (2|1), (51) 
and is expressed as a sum of two functions in the two superconductors; 

^ = ^ sl |l) + ^ 2 |2) (52) 
15 



•'max v^iingj 




FIG. 5: Sketch of a SQUID structure and a plot of the maximum current through the SQUID 
structure Jma^ ID = Max{J£^, ^1 vs the magnetic flux through the ring $ r i ng . ^junc is taken 



to be $j unc = 0. l$ ring . 



1/2 1/2 

where = p/ e~ lXl and \I/ S 2 = P2 e~ lX2 are "wave functions" for the two superconductors 
in the junction. 

By expressing the "Hamiltonian " as 

^ + g^ = ^e^^ 2A - dr |l>{2|+Te^^ lA - dr |2><l| + (g^ 1 + / x 1 )|l><l| + (^ 2 + > u 2 )|2><2|, (53) 

we obtain the model Schrodinger equation given in Eqs. f|T8|) and f|T9|) ; 



at 
at 



Then, the current is given by 



(54) 
(55) 



J = Jn sin 1 



(56) 



with S given by 



(57) 



By taking the same procedure as before, the AC Josephson frequency in Eq. (JTBT) with g = e 
is obtained. 
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VII. SQUID OSCILLATION IN THE SPIN- VORTEX SUPERCONUCTIVITY 



We now show that the spin- vortex superconductivity explains the oscillation of the max- 
imum current through a ring-shaped junction with respect to the external magnetic field. 
Let us consider the SQUID structure with two identical SIS junctions shown in Fig. The 
total current is a sum of currents through the two junctions in the ring. If the magnetic flux 
through each junction area is neglected, the total current is calculated as 



J = Jq sin 4> s + J sin 



vr$ ring /$o + w 



ring 71 ",) i 



(58) 



where $rmg is the flux through the ring region given by 



A ■ dr. 



Grin 



and w r i n g is an integer called the winding number of \ along C ring given by 



w 



ring 



2tt 



Vx • dr. 



(59) 



(60) 



Crfn 



The current in Eq. (15811 exhibits an interference pattern that depends on $ r ing and w r i ng . 
The phase difference is composed of the contribution from the Aharonov-Bohm phase from 
the magnetic filed — vr^rhig/^O; and that from the winding number of \ given by w r i ng vr. 
Note that in the Cooper-pair-tunneling case, the former is given by — 27r$ ring /$o> an d the 
latter is a multiple of 2n (i.e., 27rw ring )^^. 

The maximum current depends on whether w r i n g is even or odd. By following the well- 
known derivations, it is calculated as 



j: 



even 
max 



«o) 



sin 



7T*;, 



2#n 



2*n 



COS 



7T 



ring 



2$r 



(61) 



if w r i ng is even; and it is calculated as 



j: 



odd 



reven 
max 



(0) 



sin 



7T*j U n 

2*n 



2*o 



sin 



7T$ 



ring 



2$ r 



(62) 



if Wring is odd, where the effect of the flux through each junction region $j unc is included. 

If we take into account the change of UYi ng , the maximum current ^ax 110 i s the largest 
of Jmax an d <^max- I n Fig- a plot of J^ UID as a function of $ r in g is depicted. It exhibits a 
profile similar to the one observed in experiments^, showing peaks separated by $o- 
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FIG. 6: a: Plot of as a function of $ ring . It shows a "0-SQUID" pattern, b: Plot of J° d a d as 



a function of $ r ing; it shows a "-7T-SQUID" pattern. $j 



nuc— 0.5$ r ing is used. If <5j 



^ring is used, 



the "0-SQUID" pattern becomes the edge-junction pattern and "7T-SQUID" pattern becomes the 



corner-junction pattern given in Ref. 
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The corresponding maximum current for the 2e-charge case is given by 



J, 



2e 



sin 



3>o 



0Q 



cos 



7T$ 



ring 







(63) 



in this case, zero current is obtained at & V i ng /§o=l/2 + n, where n is an integer in contrast 
to the e-charge case; therefore, the absence of these zero-points may be used as a support 
for the e-charge order parameter. Actually, the absence was observed in the experiment, 
but the reason for it has been attributed to a background current. 

It is worth noting that if the even-w riri g and odd-w r i ng cases are separately realized, so- 
called, the "0-SQUID" and "7T-SQUID" patterns seen in cuprates are obtained, respectively 
(see Fig. [6])^i. The appearance of these patters are currently attributed to the d-wave 
symmetry of the order parameter; but it is explained due to the change of w T i ng in the 
present theory. 
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VIII. A POSSIBLE RELATIONSHIP BETWEEN THE SPIN- VORTEX SUPER- 
CONDUCTIVITY AND THE BCS THEORY OF SUPERCONDUCTIVITY 



The present re-derivation of the AC Josephson frequency indicates that the BCS theory 
based on pairing electrons is incomplete as to the origin of the persistent current; persistent 
current carriers are electrons, not Cooper-pairs. Here, we try to relate the spin-vortex 
superconductivity to the BCS superconductivity by reinterpreting the phase of the order 
parameter; we assume that Cooper pairs are broken in the surface region and spin-vortices 
are formed. 

Let us consider the BCS state vector that for the situation where the phase 9 varies over 
the sample as 

I^bcs) = J] ( Mak + ei% ak4 T ci_ ki ) |vac). (64) 

a,k 

Here, the whole system is devided into small cells, and each cell (indexed by a) is assumed 
to be described by the BCS-type state vector given in Eq. ff25l) 4 . The real parameters u a k 
and t> a k satisfy u 2 ^ + v 2 k = 1. 

In the conventional theory, 9 is thought to arise from the fluctuation of the number of 
Cooper pairs. The equation to obtain 9 and p is the Ginzburg-Landau equation^ that are 
derived from the BCS theory^; however, this derivation yields charge 2e for the carrier, 
thus, contradicts the observation that the charge on the carrier is e. In the following, we 
derive the equations of motion for 9 and p in a different maner, and relate 9 to \ that arises 
from the spin-vortex formation. 

We employ the following Lagrangian 

d 

^BCS = (*BCS|^^ - #BCS - e ^ ^a4k CT C ak<T|*BCs) 

a,k,<r 

= E - n BC s ~ 2e E Pr>a, (65) 

a a 

where pP air = X^k M ak ^ s ^ ne numDer °f Cooper pairs in the cell a, and Hbcs — 
(^bcsI-^bcsI^bcs) is the energy functional for the BCS Hamiltonian. From this Lagrangain 
the following equations of motion are obtained; 

• 2e 1 5H BC s (Rrf . 

9 + —ip = --— — — . (67) 
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We assume that even some of Cooper pairs are broken and spin-vortices are created, we 
may still use the BCS-type state vector as an approximate one if the phase % is included. 
This means that an approximate state vector is written as 

I^BCSM) = IJ («ak + ^akaik&lk) |vac), (68) 
a,k 

where new operators a a k and 6 a k are defined by 

(69) 

2 J \c a -ki J 

which corresponds to the transformation in Eq. (J!T) 
The above state vector is actually rewritten as 

I^bcsM) = J] ( Mak + e " iXa ^k4 T ci-k|) I v ac). (70) 

a,k 

The comparison with Eq. gives 6 = —\. By employing this reinterpretation many 
practical problems may be handled. However, the true equations of motion are given in 
Eqs. (HO]) and (dU), not in Eqs. (|66|) and (ISTj). 

It is known that the transition temperature T c calculated as the temperature for the 
Cooper-pair formation explains T c in many superconducting materials. However, the present 
work indicates that the Cooper-pair is not the persistent current carrier. In this respect, it 
may be worth noting that a study of correlations between normal-state properties and the 
occurrence of superconductivity indicates that properties assumed to be important within 
BCS theory rank lowest in predictive power regarding whether a material is or is not a 
superconductor.— 

The above facts may be explained if we assume that fragile germs of superconductiv- 
ity already exist above T c , and whether they exist or not is related to the occurrence of 
superconductivity. In this case, the transition temperature is the temperature where the 
germs become stable and coherently connected over the sample, that presumably occurs 
when the energy gap is formed. In this respect it is worth mentioning that the germs of 
superconductivity (or spin-vortices) are speculated to exist in the pseudogap phase above T c 
in cuprates 9 . The spin- vortices may be stable only in the surface region in the conventional 
superconductors, but even in the bulk in the cuprates due to the two-dimensionality of the 
conduction Cu0 2 plane. The anomalous pseudogap phase and the high T c value in the 
cuprates may be related to the existence of rather stable spin-vortices in the bulk 9 -. 
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IX. CONCLUDING REMARKS 

We have shown that the AC Josephson frequency 2\e\V/h actually means that the charge 
on the superconducting order parameter is e. The conventional interpretation that it is 
the evidence of the tunneling of 2e-carriers is reversed. The present result indicates that 
although the BCS theory based on the Cooper-pair formation explains bulk properties of 
superconductors very well, it fails to explain the persistent current flow which occurs pre- 
dominantly in the surface region. We speculate that some of the Cooper-pairs are broken 
or Cooper-pairs are not formed in the surface region, and spin-vortices are created there; 
then, the phase of the order parameter arises as a Berry phase effect due to the spin- vortex 
formation. This also indicates that Cooper pairs are not an absolutely necessary ingredient 
of superconductivity 9 . 

We would like to mention that there are some experimental and theoretical indications 
that spin-vortices exist in high-transition temperature cuprate superconductors 9 -. The high 
T c observed in them may be related to the fact that the two-dimensionality of the current 
flowing Cu02 plane is advantageous for the stabilization of spin-vortces. 

In the superfluidity of Bose particles, equations of motion in Eqs. §E§ and are used. 
In this case, however, the phase of the order parameter cannot be attributed to spin- vortices 
since the spin degrees-of- freedom is absent. Nevertheless, such equations are obtained by 
using two components that arise in this system due to the Bose-Einsein condensation. The 
one of the two components is the condensate state and the other is the non-condensate 
states; the phase of the order parameter may be attributed to the relative phase between 
them. We will briefly describe this possibility below. 

It is known that the ground state vector for the dilute Bose gas is given by 



where a is the creation operator for the condensate state, and describes that for other 
states; Ck are parameters and N is the total number of the particles 25 . 

We divide the whole system into small cells, and each cell (indexed by a) is assumed to 
be described by the above state vector. When the phase of the condensate wave function 




(71) 
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varies relative to the rest, slowly, over the cells, we have 

l*iv[<9]> = const. J] \al al e i2d * - ^ c^o^aO |vac), (72) 

where 9 a and iV a are the relative phase and number of particles in the cell a. 
Let us employ the following Lagrangian, 

C B = (* N [6]\ (iU^ - H^j \^n[0}), (73) 

where Hg is the Hamiltonian for the Bose system. We denote the energy functional of the 

system as 

HB = (y N [6]\H B \y N [0]), (74) 

and the number of particles in the condensate state in the cell a as 

p a = (^ N {e]\al a a0 \^ N [6]). (75) 

From the Lagrangian in Eq. (1731) . the following equations of motion are obtained; 

. 16Kb 1_ 5H b 

p = %-w = -* v ' Jw> (76) 

The equations of motion in Eqs. (176)) and (177)1 correspond to those given in Eqs. (JHJ) 
and Q). The equation ([75]) indicates that superfluid velocity is proportional to V#; if we 
neglect the non-condensate component of the state vector in Eq. (1721 . the phase 9 is the 
phase of the condensate wave function. Then, the present superfluid velocity agrees with 
that is given by A. J. Leggett 26 . In the present view, the two superfluid phenomena are 
unified in the common mathematical structure where the particle states are locally given 
in two components; spin degrees-of-freedom for the superconductivity case, and condensate 
and non-condenstate states in the Bose-Einstein condensation case. This view is markedly 
different from the view that both are due to the U(l) gauge symmetry breaking^. 

In conclusion, we have derived the AC Josephson frequency in a gauge invariant manner 
with including a battery contact effect. The result indicates that the persistent current is 
a flow of electrons, thus, reverses the current understanding of the origin of the phase of 
the order-parameter. We have presented an alternative origin of the phase of the order 
parameter; we have argued that the Berry phase that arises from spin- vortices is the origin. 
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